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Truly stable metastable states are an artifact of the mean-field approximation or the zero tem- 
perature limit. If such appealing concepts in glass theory as configurational entropy are to have 
a meaning beyond these approximations, one needs to cast them in a form involving states with 
finite lifetimes. Starting from elementary examples and using results of Gaveau and Schulman, we 
propose a simple expression for the configurational entropy and revisit the question of taking flat 
averages over metastable states. The construction is applicable to finite dimensional systems, and 
we explicitly show that for simple mean-field glass models it recovers, justifies and generalises the 
known results. The calculation emphasises the appearance of new dynamical order parameters. 



I. INTRODUCTION 



Slowly relaxing systems such as glasses or compacting granular media can be viewed as having fast, local, quasi- 
equilibrium dynamics, plus a slow, nonequilibrium drift. These two superposed motions can take different forms: 
'cage' vibrations plus structural rearrangements in glasses; bulk fluctuations plus domain-wall motion in coarsening 
problems; etc. At a given long time, the fast motion covers a region of phase space which one may picture as a 
'metastable state'. 

Even though metastable states are a familiar and appealing concept, it turns out that defining them in an unam- 
biguous manner in non-mean field models is quite subtle. This in turn has as a consequence that such standard ideas 
in glass theory as 'configurational entropy' (related to the number of metastable states) are not only hard to calculate, 
but are indeed, with the exception of some fortunate cases, approximate as concepts. 

In this paper we show how these questions can be put on a well defined basis using a formalism |jj that does not 
rely on specifically mean-field concepts. 

First of all: given that one can simulate and in certain cases calculate analytically the complete history of a sample 
starting from a quench, why should one have any need to introduce the apparently unneeded notion of 'metastable 
state' ? Indeed, these states only come into play when one wishes to make arguments such as: 'phase space contains 
such and such a distribution of states, which will be accessed with such and such a probability by a typical dynamical 
history. Long-time out of equilibrium observables can be directly calculated by averaging the observables over some 
subset of states — and further reference to dynamics may be omitted '. 

This kind of 'ergodic' argument was pioneered by Edwards who proposed that in compacting or slowly flowing 
granular systems one can obtain the correct dynamical observables by averaging the values they take over all blocked 
configurations of a certain volume. It later turned out that mean-field glass models relaxing at zero temperature 
had exactly Edwards' ergodicity property j^]: at long times any nonequilibrium observable is correctly given by the 
typical value it takes over all local energy minima of the appropriate energy density. 

A first problem arises when one wishes to apply this concept at finite temperatures (or vibration, in the case of 
granular media). There again, the mean-field case offers a suggestion: at non zero temperature Edwards' argument 
works as well, provided one substitutes 'energy minima' by 'free-energy minima' ('states'). This construction is 
possible because within mean-field we have a well defined notion of free-energy landscape, whose local minima are in 
some (but not all) cases related to completely stable distributions. However, as discussed by Franz and Virasoro [||, 
one needs to consider 'quasi-states' with finite lifetimes in order to understand the situation at finite waiting times. 
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In finite-dimensional problems, a high-lying metastable state cannot have an infinite lifetime: there is always a 
finite probability of escape through the nucleation of droplet of a more favourable phase. Hence, which distribution 
one considers as metastable depends always on which lifetimes one is considering. For example, the concept of 
'configurational entropy' (the logarithm of the number of states), ubiquitous in glass theory, has in finite dimensions 
only a meaning with a timescale attached. Moreover, even the mere definition of an Edwards distribution is not as 
simple, quite apart from the question of the validity of the ergodicity-like hypothesis it assumes. 

In section 2 we shall review the notion of metastable state within mean-field glass models and how the knowledge 
of their distribution allows in certain cases to reproduce some results obtained from the full solution of the out of 
equilibrium dynamics. We shall also mention some limitations found even at this level of the identification 'free-energy 
minimum ~ stable state'. 

In section 3 we discuss a strategy valid in any dimension for the definition and calculation of metastable states based 
on the evolution operator developed by Gaveau and Schulman We discuss how one can thus recast the question 
of configurational entropy and Edwards' distribution in a form relevant in finite dimensions at nonzero temperature 
(or stronger vibration, in the case of granular media) by considering finite-lifetime metastable states, in the spirit of 
the 'quasi-states' discussed in [||. 

In section 4 we apply this method to a simple mean-field glass model. We show how one can rederive in this way 
both the number of states and the dynamics inside a state, within a framework whose applicability goes beyond 
me an- field. 



II. A FORTUNATE CASE: MEAN-FIELD MODELS 



Consider the mean-field model of fcrromagnct: 



27V- 

*,3 



where the sum is over all spins. The spins can be Ising Si = ±1 or spherical ^ Sf = N. One can easily obtain a free 
energy in terms of local magnetisations mj = (Si) (where (•) means the average over the Gibbs measure): 

/({"*»}) = -l^^2 m * m j - h^2m l -TS{{m l }) (2) 

ijtj i 

where T = denotes the temperature and S({rrii}) is the usual entropic term: 

S(W) = ~Yj2 [(1 - mj) k(l - fflj) + (1 + nn) h(l + mt)] for S t = ±1 

i 

= — ln(l — q) spherical model (3) 

and Nq — J2iLi m f • 

The states are represented by the minima of /. At T < T c there are two, and the deeper one dominates the 
Boltzmann distribution. If h > one of the states becomes metastable: within mean-field its lifetime is infinite, but 
in finite dimensions it will decay. 



The TAP approach 



Next, consider the glassy Hamiltonians: 

E = -~i E Jn,.., p Sn-S* P (4) 
ii,...,ip 

where the Ji lt ,., t i p are independent random variables with variance p\/2N p ~ 1 . The Ising version with p = 2 is the 
Sherrington-Kirkpatrick model, the mean field version of spin-glasses. The models with p > 3 are instead systems 
having the behaviour resembling structural glasses MM ■ 
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It turns out that one can find for these models a free energy function analogous to (g) (the 'TAP' free energy M), 
in terms of local magnetisations m, : 

fTApdm,}) = -1 £ J it , ... iipmii ... mip -^[l- g f- OT f-i(l- g )]-TS({mO) (5) 

ii,...,i p 

where iVq = m? and ^({mi}) is given in (||) for the spherical and Ising cases, respectively. 

Within the TAP approach, one signature of the glass transition is the fact that the free energy (||) has, below a 
critical temperature, many (~ e aN ) minima 0,01130. The main difference between mean-field versions of spin 
and structural glasses is seen, in the TAP approach, in the way states are separated. 

Given the analogies with the ferromagnet, it seems very tempting to attribute to the minima of (^) a dynamical 
meaning of 'state'. For the models with p > 3 this has been done [ fujjl5fl by starting from a configuration where the 
coordinates are as close as possible to having a given mf, and then checking that the subsequent dynamical evolution 
is a stable, quasi-equilibrium situation confined to a region of phase space in such a way that < Si >time-average = m t- 

On the contrary, for models like the Sherrington-Kirkpatrick model the identification of all TAP-minima with stable 
states seems to breaks down. Let us formulate an heuristic argument to see this. Decreasing temperature, minima 
split in a second-order transition manner. At least a fraction of the minima are 'born' this way ]l6| ], and to get an 
exponential number of minima one needs that on average there is a division every 0(1/ N) change in temperature. 
Hence, a fraction of TAP solutions are just 0(1/ N) below their critical temperature, and under those circumstances 
barriers cannot be large enough to dynamically separate them from their 'twins'. 

FDT-temperature. Long-time out of equilibrium observables 

The dynamics of model (^) following a quench below the critical temperature can be solved analytically Jl7],[iJ . One 
finds that the system never equilibrates, and remains aging just above a threshold level of energy density eth and 
free energy density fth higher than the equilibrium ones. Given any two observables A and B, one can define their 
dynamic correlation function: 

C AB (t,t') = (A(t)B(t')) (6) 
and the integrated response xab to a field Kb conjugate to B acting between times t 1 and t: 

5(A(t)) 

XAB{t,t) = — 7 

dn B 

where the averages are over the dynamical realisation. 
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FIG. 1. A fluctuation-dissipation plot. The straight line to the left defines the effective temperature. 
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If one makes a parametric plot of XAB{t, t') versus Cab (tj t') one obtains for model (Q) with p > 2 at long times a 
curve like Fig. 2. One has in addition to the straight line with gradient — 1/T (as in equilibrium), another straight 
line of gradient, say, — 1/T e f f, associated to the slow relaxation. T e ff is the same for every pair of observables A and 
B, and it can be shown that it satisfies all the properties of a true temperature. 

The appearance of a temperature in a system is an indication of some form of ergodicity, in this case clearly not 
the usual Gibbs-Boltzmann equilibrium at temperature T. Indeed, soon after the dynamical solution was obtained, 
Monasson and Virasoro Jl^,^] observed that the temperature T e // could be reobtained from the TAP approach 
without making any reference to the dynamics: Defining the complexity (or configurational entropy) S(f) as the 
logarithm of the number of TAP solutions of a given free energy, one checks that 



dS(f) 



Teff df 



(8) 

f—fth 



How this equality follows from Edwards' assumption is discussed in 

Furthermore, one can also see that the long-time values of macroscopic observables are given by the flat average of 
their values taken over all TAP solutions of the dynamical energy density eth- Hence, we have a strong indication for 
a measure a la Edwards, this time applied to TAP states. 

If the system has large, but finite size, it will slowly approach equilibrium. In that case, a plot like Fig. 2 will show 
that the two tracts slowly tend to become parallel and T e ff(t) (now a function of time) tends to T. Inspired by work 
of Bonilla et. al. ^l), Nieuwenhuzen |2^] conjectured that one could extend the relation (||) for all times, using the 
TAP solutions of the energy level appropriate at each time. 

Later on, a two-temperature \ versus C plot (and hence the existence of an FDT temperature Teff) was seen 
to occur in realistic models such as three dimensional Lennard- Jones glasses p3|-p6[. Several of these simulations 
where performed at temperatures above the putative equilibrium glass temperature, so that the existence of a (slowly 
evolving) well-defined T e ff is surely a dynamical phenomenon, unrelated to the structure of equilibrium states. If one 
wishes to consider this as a symptom of slowly evolving flat distribution between metastable states, one finds oneself 
in the embarrassing situation that it is not entirely obvious what one means by 'metastable state' in finite dimensions, 
as we now discuss. 



The problem with finite dimensional and driven systems 



In finite dimensions and non-zero temperature nucleation arguments suggest that a distribution with dynamical 
free energy density (to be defined below) higher than the equilibrium one should decay through nucleation in finite 
times. We are hence in a situation in which we have no absolute notion of state without making reference to a 
timcscale (and hence to dynamics): two different distributions may be confused into a single state or be treated as 
two separate entities depending on whether the time to go from one to the other is smaller than or larger than the 
timescale considered. 

If we are interested in systems driven by shear or by vibration, we have the additional problem that even in the 
mean-field case the distribution is not Gibbsean within a state. In a vibrated case, the notion of stability must be 
substituted by the notion of periodicity, so that a 'state' will turn out to be a structure periodic in time. 

Before entering into the present approach, let us mention that a pragmatic way of dealing with these difficulties, 
at least at very short time-scales, is the so-called 'inherent structure' construction P7| . Though it does not solve the 
questions of principle mentioned above |28f| , it offers a practical way around applicable to concrete problems. 



III. DYNAMICAL DEFINITION OF METASTABLE STATES 



Let us consider a system evolving with stochastic dynamics, which for definiteness we shall consider is of the 
Langevin form. The probability distribution will evolve according to: 

dP ^=-HP(S,t) 



dt 
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where H is the Fokker-Planck operator. The potential energy E can be time-dependent, and furthermore one can 
add to (§§-) forces that do not derive from a potential. 

Given a probability distribution P, one can define a dynamic free energy 



F(t)= / (TP(S,t)]nP(S,t)+E(S)P(S,t)) dS 



(10) 



If H is time-independent, any stationary configuration satisfies 



HP, 



stationary 



= 



(11) 



Moreover, writing any distribution as P(x, t) = ^ Ci(t)ijji(x) where ipi(x) are the right eigenvectors of H: 



H\ij)i) 



(12) 



the evolution equation 



@ implies: 




(13) 



We see that if P is to vary slowly it has to be concentrated on eigenvectors with low eigenvalues \. Indeed, each \ 
is an inverse timescale [^9|,0 . 

In the following subsections we shall motivate and discuss an identification of the set of small eigenvalue ip^s with 
metastablc states. 



Consider first the system of two hard disks performing Langevin dynamics in a box (Fig. |2|). Clearly, if the 
disks are really not interpenetrable, there are two different ergodic components, each composed of the mirror image 
of the configurations of the other. Symmetry implies that the spectrum of the Fokker-Planck operator is doubly 
degenerate. The two lowest (zero) eigenvalues correspond to two stationary distributions. One can construct an 
associated eigenvector as the flat distribution over all pairs of coordinates of the centers of the disks such that they 
do not superpose and such that the disk 'A' is to the right and 'B' to the left, and similarly a second eigenvector 
corresponding to having the disk 'B' to the right and 'A' to the left. These are the 'pure states': any linear combination 
of these two distributions will be an 'impure' state. The next higher eigenvalues are equal to the inverse of the time 
needed for the particles to explore their ergodic component. Note, in passing, that this hard-spheres system is the 
typical example in which the inherent structure construction is not meaningful while the present one has no problem. 



In this example we have strictly two ergodic components, corresponding to twofold ground state degeneracy of the 
Fokker-Planck operator. This is indeed very general: suppose we have p ergodic components C\, ...,C Vl with typical 
times t\, t p required to explore each component. We can construct an independent eigenvector with zero eigenvalue 



Motivation 





FIG. 2. A system with two ergodic components. 
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using the stationary distribution P a (x) restricted to each component C a . These completely span the zero eigenvalue 
subspace. To show this, we calculate: 



tr[e' t ' H ] = f dx f dy(x\e~^ H \y){y\e- 1 ^ H \x) = / dx [ dy {x\e~ t ^ H \y){y\e- t i- H \x) 

= E/ dx I dy(x\e-^ H \y)(y\e-^ H \x) 
a=1 Jx£C a Jyec a 



(14) 



If we now take t* much larger than all the tj, (y\e * 2 H \x) ~ P a {y)i the equilibrium probability for y restricted to the 
crgodic component a to which x belongs. Hence: 



tr[e- 4 * H ]~^^ / dxf dyP(x)P 

a=l 6=1 J ^ C » 



(y) 



(15) 



This shows that the number of states 'below the gap' coincides with the number of ergodic components. 

In the preceding examples the ergodic components are strictly separated. However, in most applications this is not 
the case: there is in fact a passage time between components that only becomes infinite in some limit. To understand 
the construction in these cases, consider a very low temperature Langevin process occurring in asymmetric and 
symmetric double-well potentials as in Figure |. On the left of the fi gure we show the lowest levels of the spectrum 
for both cases, and at the top the corresponding eigenvectors. For the asymmetric case, the two lowest eigenvalues 
are separated by the inverse Arrhenius escape time from the highest minimum. All other eigenvalues are much higher 
(0(1)), and include the escape time from a maximum, etc. The eigenfunction labeled a is essentially positive and 
represents a "pure" state Pi, while one can make a linear combination P 2 = (a + b) that will also be positive and 
concentrated on the metastable minimum. For the symmetric well the situation is similar, but now it is the linear 
combinations Pi = (a + b) and P 2 = (a — b) that play the role of "pure" states. Any other combination of the form 
yPi + (1 — y)P2 with < y < 1 will give an 'impure' (almost) steady state. 

Note that these definitions make sense in the time-window in which we can consider the exponential Arrhenius 
times much larger that any other time involved (0(1)), and this will happen only in the low temperature limit. 

If the temperature is non-zero, a separation of timescales can happen as a result of the thermodynamical (or other) 
limit. For example, it is easy to see that a mean-field ferromagnet at < T < T c will have a similar spectrum, but 
with N playing the role of large parameter instead of l/T. The case of finite dimensional ferromagnets is slightly 
more subtle: we have there a timescale for domain excitations within a state that can be as large as a power law in 



N (the time it takes for a large domain to collapse), and a much longer timescale (~ e 
phase to the other. 



cN (d-i)/d 



for going from one 




FIG. 3. A sketch of eigenfunctions and spectrum of the Fokker-Planck operator at low temperature . 



The construction of Gaveau and Schulman 



In general, the low eigenvalue spectrum can correspond to more than two eigenvectors. One can now ask in general 
whether all these eigenvectors (or combinations of them) represent positive, stable distributions, and whether one can 
construct as many pure states as there are low eigenvalues. In a series of papers Gaveau, Schulman and Lesne 
showed this to be the case, gave recipes for the explicit construction of the pure states, proved their unicity and 
exploited this construction to study metastability. 

We shall be only descriptive here, we refer the reader to the references for the proofs, as well as other investigations 
concerning metastability. Consider for example a Fokker-Planck operator having the lowest p eigenvalues Ai, A p 
separated by a gap from the others, i.e. one can find a t* such that one can consider that: 

t*Xi « 1 for i — 1, ,..,p 

t*Xi » 1 for i = p+ 1, ... (16) 

In the previous simple example any t* such that t* is positive will satisfy this as T — > 0. The meaning of t* is clear: it 
is a timescale much longer than the relaxation into states, but much shorter than transitions between states. Clearly, 
the operator exp[— t*H] is essentially a projector onto the space 'below the gap' (up to terms of order exp[— t*A»], 
with i < p). Within the same accuracy, one can then find a basis of p right eigenvectors \Pi) which are: 

• positive: Pi(x) = (x\P; L ) > 0: 

• almost stationary: H\Pi) ^ V i = 1, ...,p 

• normalised and not zero in non-overlapping regions of space. 

The last property is related to the fact that one can also find a basis of p approximate almost stationary ((Qi\H ~ 0) 
left eigenvectors (Qi\, such that each Qi is essentially one within the support of Pi(x), zero everywhere else and satisfy 
the orthogonality and normalisation conditions: 

(Qi\Pj) ~ 5 l3 (17) 
Given any observable A, we can calculate its average within the state "i" as: 

(A)i = (Qi\A\Pi) (18) 

One can also write approximately: 

e-™ ~J2\Pi){Qi\ (19) 

i 

As a consequence the \Pi) vectors have all the good properties to represents metastable states: they are positive 
normalised distributions, non zero only on different regions of the configuration space and they are stationary on time 
scales less than t* . In the proof, as in the simple example of the previous subsection, the definition is unavoidably 
linked to a timescale: if one considers really infinite times, before any other limit, then the distinction between states 
vanishes. Furthermore, it is assumed the number of states so defined remains finite in the thermodynamic limit. We 
believe that going to situations in which this is not the case (as we will below) is indeed not entirely innocent, but 
is at the heart of quite a few problems associated to the definition of "state" in glassy systems. (We have already 
encountered such subtleties when we discussed TAP minima in the Sherrington-Kirkpatrick model) . 

Driven systems. 

The construction described above is not limited to purely relaxational or time-independent systems. Consider for 
example the case in which a system is periodically 'vibrated' or 'tapped'. One can still try to look for stationary, 
or rather, periodic situations. One can repeat essentially the same argument by considering the evolution operator 
through one cycle: 

U = T [ dt e~ tH{t) (20) 
./cycle 

where T denotes time-order. One has to look now for the eigenvectors of U whose eigenvalue is close to one. 
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Similarly, one can also work with systems driven by constant nonconservative forces (as in a sheared fluid), and 
with nonlinear space-dependent friction (as in granular systems). 

Pure barriers. 

Given that pure states can be viewed as playing the role for finite temperature that certain energy minima play for 
zero temperature, one is naturally led to ask which distributions play the role of barriers (or in general saddle points) 
in finite temperature. For Fokker-Planck processes this can be done naturally starting from the states constructed as 
above. In Appendix A we show how this can be done. 

We remark that a solvable example, in which saddle points play an important role in the spectrum of the evolution 
operators, is the Glauber evolution of the completely connected Ising model and its generalizations |30|. It has 
been found that the lowest eigenvalues, i.e. the longest relaxation times, are gathered in families, each one being in 
correspondence with a stationary, not necessary stable, point of the static mean field free energy. 



IV. FLAT DISTRIBUTION OVER STATES AT FINITE TEMPERATURE OR VIBRATION. A 
DEFINITION OF CONFIGURATION AL ENTROPY. 



Motivation. 

As we mentioned in Section II, the fact that one finds a two-temperature behaviour in mean-field glasses can be 
seen as suggesting the relevance of a measure consisting in summing over metastable states of given energy (or free 
energy) with equal weight. Armed with the construction for metastable states we discussed above, we shall see how 
this measure can be expressed in finite dimensional, or periodically driven systems. 

Expressions. 

Consider, in the spirit of Edwards' distribution, an average of an observable over states, each measured with equal 
weight. This will be relevant for the long-time out of equilibrium dynamics under the assumption that almost all 
states of given energy have the same basin of attraction. 

Using the same notation as in the previous section, we define the average over the measure ((•}} of an observable A 
in the following way: 



((A)) t .=l y £(A) i (21) 

where the subindex i* reminds us of the fact that states are now defined according to their timescale. We have: 

p p p 

= £<Qi|A|P,) =^2tr[\Pi)(Qi\A] = tr[e- ttH A] (22) 

i— 1 i— 1 i—l 



where we have used (J18|) and (19). Hence 



Note that once written this way, all reference to pure states has disappeared, except indirectly in the value chosen for 
t*. 

We often need an equation like (|23|), but restricted to states having a certain energy, particle number, etc. In that 
case we generalise ( p3| ) as, for example: 

tv[6(E - Eq) e-*' H A] t, Eo [e^ H A] 
{{A))t * {Eo) tr[J(£?-£? )e-«-*] " tr Eo [e~^] (M) 
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where tTE denotes a restricted trace. 

Once we make the assumption that all states of the same energy (or particle number, etc) have an equal weight for 
the purposes of calculating a dynamical observable, it becomes meaningful to count their number at given energy, the 
configurational entropy: 

S t *(E ) = lntr[8{E-E ) e^* 11 ] = lntr^e - ***] (25) 



If t* — > we get the microcanonical measure (and entropy, up to irrelevant constants |31|), and if we let t* — > oo we 
find no high-lying metastable states at all in finite dimensions. The dependence on t* is hence unavoidable if one is 
to obtain a finite configurational entropy in that case. 

Equation ( p5| ) defines the timescale-dependent configurational entropy. One also needs the average entropy within 
a state Sf, and the corresponding average free-energy of a state ft* = E a — Tsf. Using the construction of (p7|), ( |T§| ) 
and (|jg|), we have: 

St * {Eo) = L /E( ^ E J*'We-^\ X >) (26) 

The meaning of this equation becomes transparent in the example of the completely separated ergodic components 
of the previous section. 

Note that also the intra-state entropy is timescale-dependent. Indeed, if we set t* — ► 0, we are defining as 'states' 
the configurations themselves |n|]. On the other extreme, if t* is longer than the equilibration time, (a/je - * H \x') 
gives the Gibbs-measure, and the intra-state entropy becomes the usual entropy. In short: changing the timescale 
both changes the number of states and their nature, hence the change in configurational and intra-state entropy, 
respectively. 

Flat measures and effective temperatures 

Suppose we have a glassy system, taken through a given thermal history (an annealing protocol) to a glassy phase 
at time t, at which time its energy is E(t) = E a (and if we let the particle number or the volume change we should 
specify also them). The assumption of typicality of metastable states, is then: 

<^»'- ^ = Xt°X-"»} ~ <A)hiSt ° ry (2?) 

where the average is over several realisations of the same protocol, ending at time t with energy E(t) = E a . 

From the Langevin point of view, the left hand side of ( p^ ) corresponds to adding over all periodic trajectories 
starting from an energy E Q with period t* . Equation ( p^ ) tells us that thermal histories give the same result as a very 
particular set of trajectories, in a manner analogous as when one represents chaotic dynamical systems by using only 
the periodic orbits. 

In the zero-temperature or in the mean- field case, we can set t* = oo in ( p7| ) , and thus select the states with infinite 
lifetime. This is Edwards' prescription for granular media (T ~ 0) and the one we discussed above within mean-field. 
In finite dimensions, where metastable states eventually nucleate, we are forced to give t* a finite value. The fact that 
the choice of t* is not unique already tells us that equation ( p7j ) will be an approximation. 

The central remaining question now is what is a reasonable value for t* . Indeed, giving a value of t* determines a 
configurational entropy Sf{E ) and an intra-state free-energy: f t *(E ) = E a — Ts t *{E ). This in turn determines a 
temperature associated with the timescale as: 



dS t * 
dft* 



(28) 



(where we have eliminated the energy in favour of ft*). It seems now natural to compare the different relaxation 
times of the correlations in a given problem with t*. For example, glassy systems can often be described with two 
timescales, a fast ('/?') relaxation and a slow, (waiting time dependent if the system is aging) 'a-relaxation' t a . If t* 
is small (t* « tp), Eqn. @ gives T eff ~ T Q. If we put instead: 

t(3 « t* < t a (29) 



9 



then Tf may be different (larger) than T. This is the temperature to compare with the one governing the relation 
between correlation and response in the regime corresponding to the a-relaxation |j2| . 

There are cases in which the 'slow' relaxations happen in several timescales, becoming more and more different as 
aging proceeds pj. Then, the definition ( p8| ) immediately yields a different temperature for every widely separated 
timescale, and this would agree with what one observes from fluctuation-dissipation relations pj. 

Time dependent order parameters 

We are interested in calculating two-time correlation functions: 

' e -{e-t)H Ae -{t-t')H Be -t'H 

where r = t — t' , and J\f is the normalisation. Cyclic permutation implies: 

< AB > (r) =< BA > (t* - t) (31) 
If H is a time-dependent Fokker-Planck operator, associated with forces deriving from a potential E, we have: 

e PE H( ,-0E = (32) 

Using (B2p in (BOI) we get the time-reversal property: 



< A{t)B{t') >= itr 



e -V-T)H Ae -TH B 



-< AB > (r) 



(30) 



< AB > (r) = tr* 

where A = {tP B A&~^ and B = (e^Be"^)* 

We shall need to consider the 
system of N degrees of freedom: 



e -(f-r)H^/3E Be -0 E) i e -rH^E Ae -0E^ 



-< BA >* (r) 



(33) 



We shall need to consider the cases in which A and B are respectively Si and Si = — Let us define, for a 



C(t - t') 


= C{t) 


R(t - 1') 


= R(r) 


D(t - t') 


= D(t) 



N • 

1 . 

N ■ 

1 

N ■ 



(34) 



Because we are considering periodic trajectories, causality is violated. This means that neither D(r) nor R(t) for 
negative r need to be zero. However, it is easy to show that if (t* — r) is larger than the thermalisation time (in 
which the system is projected to its Gibbs measure), then we recover D(t) = and R(t)R(—t) = 0. The existence 
of solutions violating causality for t* large is then a symptom of large equilibration times, i.e. of glassiness. 

Using the time- reversal equation (|33|), it is easy to derive a non-causal form of FDT, valid for any t*: 



d_ 

di 



7 C(t-t')=T{R(r)-R(-r)} 



(35) 



The physical meaning of D, unlike that of C and R, is unfamiliar. If we couple the states to a time-dependent, 
random magnetic 'pinning fields' hi such that hi (t)hj (f) — F(t, t') 5ij , the fields will make a change in the number 
of metastable states, and it is easy to show that: 



D(t,t') 



1 



^ln{tr[e— ]} 



NSF(t,f) 

It is hence clear that systems with a finite number of states will have D — in the long-time limit. 



(36) 
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V. THE CALCULATION 



In the following we apply the theory discussed above to a simple glassy system: the spherical version of model (Q) . 
These models are thought to be mean-field versions of structural glasses, we shall not deal in this paper with models 
corresponding to spin-glasses for which, as mentioned above, we do not expect the present computation of states to 
be in correspondence with the TAP-equation based calculations in the literature. 

The trace of the Fokker- Planck operator, at a fixed energy density E, can be written as a functional integral |35| ] 
over the spin fields Si(t) and the response fields Si(t) with periodic boundary conditions on Si(t). Once the average 
of the trace has been performed p3, the action depends on the fields Si(t) and Si(t) through the two time functions 
C(t,t'), R(t,t'), D(t,t') only. As a consequence one can integrate out the fields Si(t) and §i(t) and get an effective 
action on two time functions: 



S/N = - [ dt (d t R(t,t') + XR(t,t') -TD(t,t'))\ t , =t+ 
Jo 



(37) 



o 



where the operator M reads: 



dt{C{t,t) - 1) + -TrlnM 
2 Jo 2 



\r — I R(t,t') C(t,t') i 



D(t,t') R{t',t) 



Since wc consider times of order one with respect to N the functional integral is dominated by a saddle point contri- 
bution. We shall obtain periodic dynamic solutions which, in the glassy phase (a) break causality, (b) have non-zero 
action, (c) satisfy time-translational invariance, and (d) satisfy time- reversal and its consequence (|35|). Note that 
(a) and (b) are properties typical of instantons, while (c) and (d) are not. In the high temperature phase there is a 
periodic solution with zero action for long times corresponding essentially to the equilibrium dynamics. 
The stationarity conditions on the action are equivalent to four equations on the two-time functions: 

C'(t) = -AC(r) + 2TR(-t) + \ [ dt"C p -\t - t")R(t' - t") + k f R(t - t")C p - 2 (t - t")C(t" - t')dt" (39) 

2 Jo Jo 

t* 



R'(t) = -XR{t) + 2TD{t) + | / di!'C v - x {t - t")D(t' - t") + k dt"C p - 2 {t - t")R(t - t")R(t" - t') + S(t) 

2 Jo Jo 

(40) 



R'(t) = -XR(t) + k f dt"D(t' - t")C p - 2 (t' - t")C{t - t") + k [ dt"C p - 2 (t' - t")R(t - t")R(t" - t') 
Jo Jo 



t 

+ k(p-2) I dt"C p ~ 3 {t' -t")R{t' -t")R{t" -t')C{t-t")- XC(t-t') + 5(T) (41) 



o 

D'{t) = -XD(t) +k f dt"D(t' - t")R(t" - t)C p - 2 (t - t") + k [ dt"D(t - t")C p - 2 (t - t")R(t" - t 1 ) 
Jo Jo 

t" 

" T>( i - i"\ T>(+" - A Kit" - +'\r.p-3(+ _ +"\ 



+ k{p - 2) / dt"R(t - t")R(t" - t)R(t" - t')C p -%t - t") - XR{t) 



n 



where k = P<,P 2 — and we make explicitly use of the time translation invariance r = t — t'. The spherical condition fix 
the value of A, which can be obtained subtracting eq. ( f4l|) to eq. ( [i"o| ) for r = 0: 

X=(p-2)(p/2 dt"C p - 1 {t-t")D{t-t") + k dt"R(t-t")R(t" -t)C p - 2 {t-t")j -2T£>(0) (42) 
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Moreover fixing the value of the energy E gives an equation on the spherical multiplier A: 

pE = -A + T(R(0 + ) + -R(CT)) (43) 

as a consequence E and A are directly related. Using the FDT relation one can show that ( f39f|42| ) reduce to a set of 
three independent equations on the functions: C(t), R(t) + R(—t) and D(t). 



A. Time-reversal, non-causal solutions 

Let us show that computing for very large t* the trace of the Fokker-Planck operator one can recover the the number 
of stable states, and the dynamics within these states. The number of stable states can be obtained by a pure static 
computation for the p-spin spherical model (p > 2) using the TAP equations [^3[ . 

For very large t* there are two possible behaviours for the two-time functions depending on the energy (and the 
model) we consider: 



if at the energy value considered there are stable states then the action evaluated in the solution has a well 
defined limit as t* — > oo. In this same limit, one expects that the two-time functions for finite r describe the 
dynamics inside a stable state (as calculated previously with other methods Jl4| ) . A careful analysis of equations 
( |39f]4"2| ) allows one to show that the asymptotic forms of two-time quantities reads for r << t* : 

C(r) = C c (r) + ^C(r), R( T )=R c {T) + ^(R(r)+r-r c ), (44) 

fl w=> + ^£) (45) 

The function R c (r) is causal, and for large r (but small with respect to t*) we have that R c (t) — > and 
Cc( r ) — > Q- Together R c and C c describe the relaxation within a state. All other functions are of order one 
when their arguments are of order one, and tend to zero when their argument is large. The corrections of order 
1/t* are not subleading in the computation of the trace, as this involves integrals over a large interval [0,t*]. 
The Edwards- Anderson parameter q and r — r c are order one constants to be determined in what follows. Note 
that the scaling of D implies that A is of order 1/t*. One can easily check that C c (r) and R c (t) satisfy the 
equations describing the equilibrium dynamics inside a state studied in ||l4| : 

C' c (t) = AC c (r) +p/2 f dt"C%- l {t - t")R c (t' - t") 



o 



+ */ Rc(t-t")C*- 2 (t-t")Cc(t" -t')dt" +p 2 q p ~ 1 (r-r c )/2 (46) 



R c (t) is causal and related to the correlation function through the FDT relation R c (t) = —1/TC' c (t)6(t). 

• if for the energy value considered there are no stable states then the behaviour of the two-time functions is 
similar to the previous one except that D has a part of order of one for finite time D Q and a part of order or 
1/t* for r - t*: 

D(t) = D (t) + t Ld(^) (47) 



(t*) Vt 

As a consequence C(r) and R(t) do not satisfy equations 'within a state', and one has to solve a set of three 
equations on C, R and D in which R is not causal also for infinite times (t*). 

In the rest of this subsection we consider an energy such that stable states exists and we compute the zero-frequency 
values of C, R and D. Using the asymptotic form introduced before one finds: 
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t* 

C{r)dT = qt* + 0(1) (48) 



f 1 

J R c (r)dr - r c = (49) 
t* ft* 

R(r)dT= R c (r)+r-r c + 0(l/t*)=r + 0(l/t*) (50) 



o Jo 

t* 



J D{T)dT=±.+0(\/{t*f) 



(51) 



Moreover for large t* the relation ( |43|) reduces to the usual one: pE 1 = — A + T/2 + 0(l/f ). Therefore we can consider 
A as a fixed parameter. 

Integrating eq. (p9) between and t* and taking the leading order in t* we obtain: 



>1 I -A + ^(1 - g)g p - 2 + ^(1 - g^ 1 ) + ^(r - r c )g^ 2 ) = (52) 



Subtracting eq. ( p9| ) evaluated in r = to ( |5^ ) we get the usual equation on A and g: 

^ =1^ + ^(1-^) (») 

Since the spherical multiplier is a parameter, this equation fixes the overlap. It is useful to write this equation in a 
way which is directly related to the static computation. Using the following notation: 

q^p£ = -A+ - q^) fcH^l _ q) z = (i _ q )gP/*-i/T (54) 

where £ corresponds to the zero-temperature or radial energy which appear in the static computation, one can rewrite 
(|53|) as the usual static equation on q |l3| 

l+pz£ +p(p- l)z 2 /2 = (55) 

From equation ( |52| ) we find the value of r: 

r = --q 1 -P' 2 £ (56) 
P 

Integrating (|4^) between and t* and using eq. ( |52] ) we get the value of <i: 

q p ~ 1 d= -2/p + 4£ 2 /p 2 (57) 

Finally we confirm, using the relationship between A and q, that the equations on C c (r) and R c (t) (|4^) are indeed 
the same ones found in [ flif for the relaxational dynamics inside the stable state with overlap q and energy E = 
(-X + T/2)/p. 

B. The configurational entropy 

To obtain the number of stable states we have to inject the solution of eqs. (|39|-^2|) into the action ( |37| ) and then 
take the long time limit. Using the compact notation Q(t; t*) = (C(r; t*), R(t; t*), D(t; t*)) for the set of the two-time 
functions, we decompose the asymptotic solution as Q(r;t*) = Qo{T;t*) + Q\(r;t*)/t* , where Q reads: 

Qo(r;t*)=(c c (T),R c (T) + ^,^L'\ (58) 

and C c (t), R c (t) are the solution of the relaxational dynamics inside a stable state (EG) with the values of r,r c ,d 
determined above. Using this decomposition and the fact that Q is a saddle point of fl37j) we find for very large t*: 
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S(Q)~S(Qo)-^®0®^ + ... (59) 



An explicit computation shows that the second term of (q9|) vanishes in the long time limit provided that the corrections 
to q are of an order less than 1 /t* . This seems natural to us since we are considering the relaxation dynamics inside a 
stable state. As a consequence the (annealed average of the) logarithm of the number of stable states coincides with 
S(Qo). When we inject Qo into ( |37|) the first two lines can be easily computed and read: 

- A(r - r c ) + \ dq?- 1 + ^{r~ r c )(l - -Z^ 1 ) + ^^(r r c )V~ 2 (60) 

whereas the computation of the third one, which reduces only to TrlnM/2, is slightly more subtle. Since the operator 
( |38| ) is diagonal in Fourier space, we get: 



Tr 1„ M : : - ln(r 2 - qd) + \ £ In ( ) , (61) 



L . 1 , , 2 n , 1 / AM C(W) 

-TrlnM = - ln(r - fld) + - > In ~ \ , 

where the Fourier transform of a function F(t) is defined as 



f lirn 
F(uj)= e- lUT F{ T )dT , oj = — n = 0,±l,±2,... (62) 
Jo t* 

The function R c (t,t') is causal and the associated operator is upper triangular with diagonal elements equal to unity. 
Its determinant hence is one (here the Ito convention is crucial) , as it should because it does not give any contribution 
to the action in the standard case. As a consequence we have: 



R(oj) C(u) 



= - ln(r c ) , (63) 



Collecting all the pieces together and using the equation on q, r, r c and d obtained in section V A we find that the 
action S(Qq) reads: 



S(Qo)/N = 1(1 + lnp - In 2) - £ 2 + \ { - -A | lu( -<f - \,'f 2 - £~ ) (04) 



Sr. — — 1 



'2(p-l) 



V p 

As expected, this expression coincides with the logarithm of the number of TAP states computed by Crisanti and 
Sommers 13 1. 

Note that this formula is correct only for £ < £ c . For £ > £ c the formalism tells us that there are no stable states 
as follows: in this energy regime D and A remain of order one even for infinite t* and the action acquires a negative 
contribution of order t* . This is as it should: since for these energy values there are only metastable states with finite 
life-times, the longer we set t* the less metastable states we find. 

In general for finite dimensional glassy systems the interesting quantity will be the logarithm of the number of 
metastable states with finite lifetime, which can be obtained plugging the solution Q(t) into the action S for a finite 
value of t* . In the following we compute this quantity at zero temperature. 



C. The zero temperature case 

At zero temperature the equation ( ficfr ) is particularly simple. Introducing the notation R(t) = R c (r) + (r — r c )/t* 
we find that the equation on the Fourier transform of R c reads: 

- iuoR c {uj) - \R c (tu) + EV-Zl ~1r c (gj)R c (uj) + 1 = (65) 
For each frequency there are two solutions: 
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Rf{u) = (-A± v/A 2 - 2p(p- 1)) , A = A + io; (66) 

In the following we focus on the two solutions R£ (r) and i?~ (r) which correspond respectively to taking the Fourier 
transform of and R~(uj). Using that at zero temperature C(t) = 1 one can decompose S(Q) in two terms 

such that all the dependence on t* is contained only in one of them: 

S/N(Q±) = 1(1 - bxp In 2) - £ 2 + 1 ( £ T ^ = ^ +ln(-f T V^^f) (67) 



p(p - 1) ''' 



jf J^=(T)^=(-r)dT + l^M^(w)4 ± (-w)) 



The computation of the second line is performed in the Appendix B. It turns out that, as in the static case |l2| ], the 
dominant contribution is given by R~ for £rsb < £ < £c, by R£ for — £ c < £ < —£rsb and for — £ c < £ < —£ c the 
two saddle point contributions are the same (see Appendix B). Note that at zero temperature, £ is the energy density 
of the system. The final result is 

S(Q) = Re f 1(1 - \np - In 2) - £ 2 + ~ ^ T ^^ —^ + H~£ T ^~£l) j (68) 
dujp p (u! +p£ ) ln(l — exp(— t |w|)) + 1 / du)p p (ui + p£)u> 

J — oo 

where = y/2p(p — 1) — x 2 /(Trp(p— 1)) is the Wigner semi-circle law. In fig. |we plot (|68| ) forp = 3 and different 

values of t* as a function of the energy density £ . For very large values of t*, S(Q) converges to the logarithm of the 
number of stable states. Note that for a finite dimensional system we expect a similar behavior but with a vanishing 
curve for infinite t* . Finally, we remark that the formula ( |S8|) has a simple interpretation. In fact the first line coicides 
with the number of saddles with energy density E. Moreover, since the spectrum of the Fokkcr-Planck operator for 
an harmonic oscillator with frequency u is Q E n = (n + 1/2) |w — uj/2, the second line of d6§| ) corresponds to the 
contribution due to a collection of harmonic oscillators with frequency distributed by the semicircle law centered in 
—p£. This distribution is exactly the same of the eigenvalues of the energy Hessian evaluated in saddles with energy 
density £ Q. As a consequence, at zero temperature, the spectrum of the Fokker-Planck operator for the p-spin 
spherical model coincides with the one obtained making an harmonic expansion around each saddle (also the instable 
ones). 




t =2000 



-0.02 1 1 1 1 1 1 1 1 1 

-1.1717 -1.1617 -1.1517 -1.1417 

e 

FIG. 4. Time dependent configurational entropy s(e, t*) at T — > for p = 3 as a function of the energy density e and of the 
time t*. From top to bottom: t* = 5, 10, 100, 1000, 2000, +oo. Note the infinite time limit curve with positive configurational 
entropy a mean-field artifact. 
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VI. THE TWO-GROUPS ANSATZ AND SUPERSYMMETRY BREAKING 



Twenty years ago, when people started to search for replica symmetry breaking solutions of the Sherrigton- 
Kirkpatrick model, Bray and Moore |55| proposed a two-groups Ansatz for the famous Q a ,b matrix j|. At a first sight 
the results seemed a little bit strange since the lim„^o Z n ^ 1 ! ( 7 means the average over disorder) . It turned out 
p6| that in the limit of n going to zero, the logarithm of Z n equals the the logarithm of the number of TAP states, 
i.e. the long-time limit of the configurational entropy (p5|). The reason of this "coincidence" has been completely 
obscure until now. For instance, for many mean-field spin glass models the complexity was computed starting from 
TAP states and, after, it was checked that the two-groups Ansatz gave back the correct result ]3^j37| , p8[ | . 
Using the properties of the dynamical solutions presented in previous sections we can unveil why the Bray and Moore 
Ansatz allows one to calculate the long-time limit of the configurational entropy (|2^). In fact this Ansatz is isomorphic 
to zero frequency part of the dynamical calculation for t* — > oo. This can be shown by the supersymmetric formalism 
for Langevin dynamics p3| , ^8[ . Within this framework all the two points correlation functions between fields Si(t) 
and Si(t) can be encoded in: 

N 

Q(l,2) = 5^(^(1)^(2)) l = (Tt 1 ,9 1 ,9 1 ) (69) 

i=l 

Si(l) = Sifa) + 1 0iS i (t) + c i (t l )9 1 +3icj(ii) (70) 



where 9±,6i are Grassmann variables and a(t),a(t) are fermion fields [ p3|j3q| . Using this formalism the dynamical 
solution giving back the configurational entropy reads at large times: 



<.H 1.2i = f,)..! 1.2) +>j (Mi - «^2)^+9iCi¥ 2 ^2 (7:1 ! 



Qc(l, 2) = ( 1 + ^(0i - 9 2 )[9i +e 2 - (9 1 - 2 )sign(ii - t^-^j (C c (h - t 2 ) - q) (72) 

The function Q c (l, 2) is supersymmetric, whereas the last two terms in the right hand side of (^Tj) break supersym- 
metry. 

On the other hand, the two-groups Ansatz consists in a symmetric Q a ,b matrix: 

B B 

Qab = lH a — b < ra, Q Q b = 1 m < a — b < n (73) 

m m 

Q ab = A-\ a^b; a,b<m (74) 

m 

Q a b — A a 7^ b; ra < a,b < n (75) 

m 

C 

Q a b — A 77 a < m, m < b < n (76) 

m l 

where one has to take the m — * oo and n — > limits. The functional dependence of the dynamical and replica free 
energy of, respectively, Q(l,2) and Q Qi h is the same Q. Indeed, the kinetic term in the dynamical free energy, 
which does not have a correspondent in the static case, is zero for the dynamical solution ( |7l| ) . Moreover if one puts 
B = (r — r c )T, 2C — T 2 d, A — q the two matrices Q(l, 2) and Q a ,b lead to the same results under tracing, convolution 
and term by term product. For instance, one can easily obtain: 

5] f(Qa.a) = f dlf(Q(l, 1)) - 2B (77) 
E = / dld2f(Q(l, 2)) = 2(/'(l) - f{A))B + f"(A)B 2 + 2f'(A)C (78) 

a, b 

Tr a , h log Q = Th, a log Q = -2 log(l - A) + log((l - A + B) 2 - 2AC) (79) 

As a consequence the computation, which make use of the two-groups Ansatz, is isomorphic to the dynamical one for 
t — > t*. Therefore the replica symmetry breaking scheme encoded in this Ansatz can be finally understood: it is a 
way to implement the dynamical computation in a replica formalism. There are, however, two important differences 
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between the two approaches. First of all, in the dynamical computation we are not free to choose between different 
Ansatze the one which gives back the long-time limit of the configurational entropy (|2^) , but we have simply to solve 
the equations of motion. This clearly makes the procedure inambiguous, unlike the case of the replica computation. 
Moreover the two approaches lead to the same results only if a dynamical solution with the correct values of q, d, r 
exists. It could then happen (cf. the discussion on the configurational entropy of the SK model) that the equations 
on q, d, r admit a solution but there is no dynamical solution corresponding to these values. As a consequence even 
if the static computation, i.e. the sum over all TAP solutions or the computation by the two-groups scheme, predicts 
the existence of an exponential number of stable states the more correct dynamical calculation does not. 



VII. CONCLUSIONS 



In this work we have shown how to put the questions related to metastable states in glasses in a manner valid for 
finite-dimensional systems. 

We have used the construction of Gaveau and Schulman to define the metastable states. This construction requires 
the existence of a 'gap' in the lifetime, so that one can associate 'states' with distributions that are stable for much 
longer than a given time t* , and 'transient processes' those that decay much faster. There is no such gap in real 
glasses, so our use of this construction has to be considered partly as a definition inspired in the cases where there is. 

A reasonable criterion for the relevance of any quantity will hence be that they are not too sensitive to the exact 
value chosen for t* . For example, if we consider a temperature associated with the a-relaxation as T t » with t* ~ t a> 
this definition is meaningful to the extent that it is stable with respect to a change in t* of, say, an order of magnitude. 

Next, there is the question as to whether T t * indeed reproduces the fluctuation-dissipation temperature. This and 
other results depend on the validity of the flatness hypotheses a la Edwards (for which positive evidence begins to 
appear J39|). In this paper we have formulated the hypothesis in a manner applicable to positive temperatures and 
finite dimensions (as well as to vibrated systems) — but we have not attempted to prove it. It may be, however, that 
writing it in the form (p7j), can be a good starting point for doing this. Moreover, the form ( p7j ) (and (|28|)) lends 
itself naturally to a generalisation to cases in which a system has more than two widely separated timescales and 
temperatures. 

Finally, the computation in section V has allowed us to check the mean-field results without relying on the TAP 
states, themselves an intrinsically mean-field concept. The kind of solutions that dominate break causality and 
have positive action, but satisfy time-reversal and a non-causal form of FDT. Unlike the barrier-crossing solution of 
Lopatin and Ioffe |io|| , they have in this sense only some of the properties of true instantons. Moreover, the dynamical 
computation unveils the meaning of the two-groups Ansatz |5(| , which allows one to compute the number of stable 
states within a replica formalism. 
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VIII. APPENDIX A. PURE BARRIERS. 



Let us show how to define 'pure barriers' using a supersymmetric extention of the Fokker-Planck operator. 
For a system with N degrees of freedom, introduce the N fermion creation and annihilation operators 0.$ and a\ , 
with anticommutation relations [aj,aj] + = Sij. Define the supersymmetric charges as: 

Q = (Tpi - iE.i)a\ ; Q = piOn (80) 

The supersymmetric operator Hsusy 

Hsusy = [Q, Q]+ = H + -J^a] ai (81) 

commutes with the charges, and with the fermion number operator. In the zero-fermion subspace it coincides with 
the Fokker-Planck operator. 

Applying the operator Q to all but the lowest Fokker-Planck (zero-fermion) eigenvectors, one obtains a one-fermion 
eigenvector. Going back to the example of low-temperature dynamics in a potential of section III, one can see that 
the lowest one-fermion eigenvectors correspond to distributions associated to the barriers. Indeed, one can convince 
oneself that in a low temperature multi-dimensional system, the lowest eigenvectors with k fermions correspond to 
barriers with k unstable directions (see Ref. where these questions are discussed in detail, and used to derive 
Morse theory results). 

Now, by analogy with the argument motivating the construction of pure states, it is reasonable to define 'barrier 
distributions' in general as the lowest eigenvectors of Hsusy with lowest eigenvalues in subspaces with k fermions. 
This definition will make sense whenever there is a gap in the spectrum, whatever the origin of such a gap (large N, 
low T, etc). One can then use traces of e~* HsuSY to calculate expectation values. 
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IX. APPENDIX B. 



The aim of this section is to calculate the last line of (^). Instead of making the computation by brute force, 
we will use the exact results that can be obtained in the p = 2 case. For p = 2 the spherical model is a simple 
collection of harmonic oscillators with frequency distributed with a semicircle law centered around the value of the 
spherical multiplier. The spectrum of the Fokker-Planck operator for an harmonic oscillator with frequency u> is [[33): 
E n = (n + l/2)|w| — u>/2. Therefore for p = 2 and at zero temperature, the logarithm of the trace of the evolution 
operator reads: 

lntTre-'*^) = -N f dujp{uj + 2S) ln(l - cxp(-i* M)) + N ( du;p(uj + 2£)t*u (82) 

J J —oo 

where p(x) is the Wigner semicircle law. One can obtain this result also by the functional formalism. In the 
computation of section ( |V CJ ) we obtained an action which for p = 2 reads: 

S/N = ±-£ 2 + ±(£ Tv / ^ _ T) 2 +In(-£TV / £ J_ T) (83) 

n n 

where we have written R(t) instead R c (t) since r = r c for p = 2. This is directly related to the fact that there is 
no configurational entropy for p — 2 and therefore the two-time functions relax asymptotically faster than 1 jt* . For 
£ < — 1 the p = 2 model is a collection of stable harmonic oscillators. One can easily compute R c (t) starting from 
the result for a single oscillator and integrating over the Wigner distribution. As expected, this function coincides 
with R~(t). As a consequence for £ < — 1 we expect that lnTre~** ff = NS{Q~). On the other hand for £ > 1 the 
p = 2 model is a collection of unstable harmonic oscillators which can mapped to the previous case changing the sign 
of each ui. As a consequence for £ > 1 we expect that lnTre~** ff = NS(Q + ). In the intermediate energy regime a 
priori one has to consider both solutions. Since the functional computation should give back the result (|8^), the last 
line of ([83] ) reads: 

~\ E R± W R± W + \ E H&W&i-n)) = (84) 

n n 

= -J dujp{uj + 2£) (-**! + In |e wt * /2 - e^'*/ 2 |) ± in J dujp(u + 2£) 
if the determination of the logarithm in the first line of ( |83|) is such that In— 1 = — in. Note that for £ < — 1 all the 



oscillators are stable, as a consequence the last term in ( p4[) vanishes and (83) coincides with (|32|). For £ > 1 all 



the oscillators are unstable, as a consequence the last term in (|84|) equals Z7r and cancels the —in coming from the 
first logarithm in (83), and one obtains the same results that for £ < — 1 with an additional term J duip(uj + 2£)t*u>. 
Finally, in the intermediate energy regime the first line in ( ]83| ) is complex and its imaginary part cancels exactly the 
imaginary contribution coming from the last term in ( |8~I| ) . In this case the the two saddle point contributions are the 
same, therefore to obtain ( |S2| ) one can consider only one of them. However to obtain the expected value of R(t) one 
has to sum on saddle points. 

For p greater than two, the equation on R c has the same form of (^0|) but for a p = 2 spherical model at zero 
temperature with a variance of the couplings J 2 = p(p — l)/2. As a consequence the last line of ( |67| ) reads: 

-J dijp((j+p£) (-**!+ In \e ut */ 2 -e- ut */ 2 \)±in J dup(u+p£) (85) 
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